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Abstract. We study the effect of thermal fluctuations on a probe qubit interacting 
with a Bose-Einstein condensed (BEC) reservoir. The zero-temperature case was 
studied in [Haikka P et al 2011 Phys. Rev. A 84 031602], where we proposed a 
method to probe the effects of dimensionality and scattering length of a BEC based 
on its behavior as an environment. Here we show that the sensitivity of the probe 
qubit is remarkably robust against thermal noise. We give an intuitive explanation for 
the thermal resilience, showing that it is due to the unique choice of the probe qubit 
architecture of our model. 
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1. Introduction 

Over the past three decades quantum computing has been the holy grail of quantum 
information sciences [1]. Lately, there has been a notable shift of focus from studies of 
circuit based quantum computers, where a long computation involving many qubits is 
broken down to elementary one- and two-qubit quantum gates, to studies of quantum 
simulators, where a physical system is modelled using another physical realisation of 
the original Hamiltonian [2]. One celebrated example of the latter is the simulation 
of a Bose-Hubbard model in optically trapped ultracold gases. The mapping between 
the Bose-Hubbard Hamiltonian and the Hamiltonian describing ultracold atoms in an 
optical lattice was proposed in 1998 by D. Jaksch et at [3] and realized experimentally 
a few years later by the group of I. Bloch [4]. Since this milestone there has been an 
explosion in studies of systems that can be simulated with ultracold quantum gases 
[5, 6]. 

Quantum simulations have been considered in the context of open quantum systems 
with proposals of simulating the spin-boson model using, for example, a quantum dot 
coupled to a Luttinger liquid [7] or a more general Bose-Einstein condensed (BEC) 
reservoir [8] . Both cases realise the independent-boson Hamiltonian with an Ohmic-like 
spectrum of the reservoir [9] . Another proposal in this direction was presented in Ref . 
[10], where an impurity atom in a double well potential is immersed in a BEC reservoir, 
forming a spin-boson model with a reservoir spectral function that can be tuned from 
sub-Ohmic to Ohmic to super-Ohmic. With a super-Ohmic spectrum the spin-boson 
model can acquire non-Markovian properties [11], hence simulating a prototype of a 
non-Markovian open quantum system model. 

Non-Markovian systems have been the subject of intense studies in recent few 
years, boosted by the recent introductions of several non-Markovianity measures that 
define and quantify the amount of non-Markovianity in a quantum process [12, 13, 14]. 
Fundamental interest in non-Markovian processes stems from the fact that Markovian 
dynamics is typically only an approximation, which is no longer valid when considering 
shorter time-scales and/or stronger system-environment couplings. Furthermore, there 
have been proposals for using non-Markovianity as a resource in the context of quantum 
metrology [15] and quantum key distribution [16], to name a couple of examples. 

Spin systems coupled to ultracold gases are not only important for quantum 
simulations, but also because they can be used to probe the ultracold gas: the 
way a spin decoheres under the action of the ultracold gas may depend crucially on 
certain properties of the gas. Hence it is possible to recover information about the 
large and generally inaccessible environment by looking at the spin only. Indeed, the 
aforementioned independent boson models can be used to probe the Luttinger liquid 
parameter [7] and the density fluctuations of the BEC [8]. In Ref. [11] we demonstrated 
that the non-Markovian properties of the impurity atom in a double well potential give 
indications of the effective dimensionality of the BEC reservoir. In this work we further 
consider this model, taking a step towards a more realistic scenario by considering the 
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effect of thermal fluctuations on the sensitivity of the probe qubit. We demonstrate 
that the double well qubit is remarkably robust against thermal noise and therefore a 
good candidate for probing ultracold gases. 

2. The model 

We consider a qubit model based on a single atomic impurity trapped in a double well 
potential, where the pseudo-spin states are represented by presence of the impurity atom 
in the left or the right well of the double well potential. The qubit is immersed in a 
thermally equilibrated Bose-Einstein condensate reservoir. The Hamiltonian of the total 
closed system is (h = 1) 

H = H A + H B + H AB , (1) 



where 



H A = J rfx 



2m A 



*(*) (2) 



is the Hamiltonian of the impurity atom with ^(x) the impurity field operator and 
Va(x) the double well potential formed by an optical lattice, 



H R = / dx $ 



^ + y B (x) + ^i>t(x)<l(x) 



$(x) (3) 



is the Hamiltonian for the BEC with $(x) the condensate field operator, Vb(x) the 
harmonic trapping potential and gs = 4:irh 2 aB/mB the boson-boson coupling constant, 
and finally 

Hab = 9ab J rfx ^ t (x)$ t (x)<l>(x)^(x) (4) 

is the interaction Hamiltonian with g^B — ^h 2 aAB/ m AB the coupling between the 
impurity atom and the condensate gas. Masses of the impurity atom and the background 
bosons are m A /B, ^ab = (ttia+tub) / (mAiriB) is their reduced mass and clb/ab are the s- 
wave scattering lengths for the boson-boson and impurity-boson collisions, respectively. 

We would like to stress that the boson-boson scattering length can be manipulated 
by Feshbach resonances, providing a controllable environment of interacting bosons. 
The significance of this is twofold: for one, many typical models of open qubit systems 
assume a non-interacting bosonic environment and it is fundamentally interesting to 
study interacting models. Secondly, the ability to have experimentally feasible and 
precise control over the environment is vital for reservoir engineering. 

We assume that the BEC is trapped in such shallow potential that it may be 
considered to be homogenous, while the double well trap for the impurity atom is so 
deep that tunneling from one well to the other is suppressed. The condensate is treated 
in the Bogoliubov approximation, assuming weak to moderate boson-boson coupling. 
After imposing these assumptions on Hamiltonians (2)-(4) the qubit dynamics can be 
derived without any further approximation. The result is purely dephasing dynamics 
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of the qubit with constant populations and off diagonal elements of the qubit density 
matrix decaying as 

|poi(t)|=e- r «|p i(0)|. (5) 
The decoherence factor is 

r (i) = 8g AB n 2^(K|-KI) e 7 ^ coth sm ( k - L )>( 6 ) 

where no is the condensate density, |ttfc| and \vk\ are the k-th Bogoliubov mode 
amplitudes with energy E\. = y2eknog~B + e|, free modes have energy = h 2 k 2 /(2iiib), 
t is the width of the impurity wavefunction, assumed Gaussian, in each well of the 
double well and /3 = 1/ksT. The spatial separation between the two wells is L. For a 
detailed derivation of the decoherence factor see [10, 11]. 



3. Non-Markovianity Measure 

With the decoherence factor at hand one has a full description of the qubit dynamics. 
We proved in Ref. [11] that in this case non-Markovianity is directly connected to 
the negativity of the decay rate 7(t) = dT{t)/dt. In this Section the connection is 
briefly reviewed and extended to the case of thermal reservoirs. We focus on the 
approach of Ref. [12], which defines Markovianity to be a property of a dynamical 
map p(0) i — y p{t) = $(£, 0)p(0) that monotonically decreases the distinguishability 
D[pi, P2] = \\pi — P2I of any two system states pi^if)- Non-Markovianity is then the 
ability of a dynamical map to temporarily increase the distinguishability of two states. 
The temporal change in the distinguishability a = dD[pi, p2,t]/dt can be associated to 
information flowing from the system to its environment (cr < 0) or back to the system 
(cr > 0). The amount of non-Markovianity in a quantum process is given by the cumulant 
of the positive information flux, M = max pi P2 J a<0 dsa(s), with a maximization done 
over all possible pairs of states to find the largest amount of information that the system 
can recover from the environment. 

The maximization required in the calculation of the non-Markovianity measure is 
generally difficult. In the case of pure qubit dephasing, however, it has been proven 
that the optimising pair is formed by two antipodal states in the equator of the Bloch 
sphere and in this case the measure can be recovered analytically. One finds easily that 
information flows back to the qubit from the environment iff the decay rate is negative. 
Moreover, in the model considered in this article there is at most a single interval of 
time a < t < b such that j(t) < and thus we introduce a normalized version of 
the non-Markovianity measure, which measures the amount of recovered information 
against the amount that was lost to the environment in the time interval < t < a. 
Summarizing, the measure we use in this work to study the non-Markovian properties 
of a qubit dephasing in a BEC environment is 

^ = e -r(0)_ e -r ( .) » 7(*) = -^<0^tGM]. (7) 
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In Ref. [11] we studied the changes in the non-Markovianity measure induced by different 
effective dimensions of the reservoir and for a range of different values of the scattering 
length of the environment, assuming a zero-T environment. We found the existence 
of a dimension-dependent critical scattering length such that when < a# < a crit the 
dynamics of the qubit is Markovian and when as > a C ru it is non-Markovian. The 
dependence on the effective dimension of the BEC is such that a crit ^D < a>crit,2D < 
0'crit,iD, that is, the higher is the dimension the smaller is the critical scattering length. 
Hence one has at hand a model where the Markovian to non-Markovian crossover can 
be controlled by changing either the scattering length of the background bosons or 
by lowering the effective dimension of the BEC. Conversely, one may deduce these 
properties of the environment by looking at the qubit only, without directly measuring 
the BEC. In this work we proceed to consider the effect of thermal fluctuations on this 
result. Thermal fluctuations can, in principle, wash out non-Markovian effects and thus 
compromise the sensitivity of the probe qubit. Fortunately we find the double well qubit 
model to have a remarkable robustness against thermal effects, as shown in the next 
section. 

4. Results 

Fig. 1 shows the decoherence factor T(t) and decay rate ^(t) for a three-dimensional and 
a one- dimensional 87 Rb condensate with a range of temperatures T = to 200 nK and 
T = to 20 nK respectively. We take the same parameters as in Ref. [11] and a fixed 
value (ib = ciRb for the scattering length. In the ID case the negative part of the decay 
rate, indicating the existence of non-Markovian effects, decreases in size with increasing 
temperature until it reaches a critical temperature of around T = 6.5 nK, where it 
vanishes completely. When this happens, the qubit dynamics is Markovian. In the 3D 
case the negative part of the decay rate splits into two lobes: the low temperature lobe, 
enclosed by the line corresponding to the zero-T decay rate, and the high temperature 
lobe, indicated by the high temperature decay rate. The transition between the two 
lobes, indicated by the decay rate for T = 50 nK, corresponds to the transition from the 
low temperature regime to the high temperature regime. We demonstrate next that the 
transition from low to high temperatures is also clearly visible in the non-Markovianity 
measure. 

The temperature dependence of the non-Markovianity measure M is shown in Fig. 2 
for all three effective dimensions. In the case of a quasi-lD condensate, the system is 
non-Markovian only for very small temperatures T ~ 1 nK and for a higher temperature 
thermal fluctuations wash out memory effects in the qubit dynamics. When the qubit 
is embedded in a quasi-2D or a 3D condensate, the dynamics is more robust against 
thermal effects. In these cases the non-Markovianity measure is almost constant for low 
temperatures when coth(/3£'fc/2) m 1. As the temperature is raised, the value of the 
non-Markovianity measure decreases as the system moves towards the high temperature 
regime, while in the high-T regime the measure increases in value again. The minima 
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Figure 1. (Color online) Top: Dccoherence factor, Y\n{t) and decay rate, 7ir>(i) 
for a one-dimensional environment, for temperature ranging between K (solid blue 
line) and 20 nK (dotted red line). The T = 6.5 nK line, shown in dashed green, 
show the transition from non-Markovian low temperature limit and the Markovian 
high temperature limit. 

Bottom: Dccoherence factor, T^£>{t) and decay rate, 7313 (t) for a three-dimensional 
environment, for temperature ranging between K (solid blue line) and 200nK (dotted 
red line). The T = 50 nK line, shown in dashed green, shows the intermediate stage 
between the low and high temperature limits. 



in Fig. 2 corresponds to the decay rate moving from the low-T lobe to the hight-T lobe 
(see Fig. 1). In the high-T regime coth(/3Ek/2) ~ (f3Ek/2)~ l , and temperature acts as 
a coefficient for the decoherence factor of Eq. (6). Consequently the whole dynamics is 
amplified, leading also to higher values of the non-Markovianity measure. 

Finally, Fig. 3 shows the non-Markovianity measure M as a function of the 
manipulated scattering length of a Bose-Einstein condensate for temperatures T = 0.5 
nK and T = 100 nK. In both cases we reproduce the main result of Ref. [11], namely 
that the qubit system has a transition from Markovian to non-Markovian dynamics 
with increasing scattering length, and that the critical scattering length depends on 
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Figure 2. (Color online) Non-Markovianity measure M as a function of temperature 
for a quasi-lD (blue solid), quasi-2D (green dashed) and 3D (red dotted) environments. 
Scattering length is fixed at as = a-m- 




Figure 3. (Color online) Non-Markovianity measure TV as a function of 
relative scattering length 0,3/ am for quasi-lD, quasi-2D and 3D environments with 
temperatures T = 0.5 nK and T = 100 nK. 

the effective dimensionality of the BEC: a crit3D < a crit ^D < a cHt,iD- For high 
temperatures the quasi-lD environment is unable to return information back to the 
system, leading to purely Markovian dynamics. However, since the quasi-2D and the 
3D environments still induce non-Markovian dynamics, the information obtained on 
the effective dimensionality of the environment by looking at the qubit dynamics is the 
same. Thus we find our main result: thermal effects do not compromise the ability of the 
probe qubit to detect information about the effective dimensionality and the scattering 
length of the environment. 
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5. Discussion and Conclusions 



The remarkable sensitivity of the probe qubit we propose in this paper derives from its 
general robustness against thermal effects. This, in turn, is due to the very specific qubit 
architecture we choose. The deep double well potential, in which the impurity atom is 
trapped, imposes limitations on the contribution of certain Bogoliubov modes to the 
qubit dynamics. The decoherence function T(t) is defined as an integral over all modes 
k, however in this model there are two cut-off momenta: 1/r, relates to the size of each 
harmonic potential in the double well and 1/L characterizes to the distance between the 
minima of the wells, and only excitations corresponding to 1/L < k < 1/r contribute 
to the dynamics. For high temperatures the temperature-dependent term coth((3Ek/2) 
diverges at Ej~ = 0, that is, when k = 0. However, the lower cut-off frequency excludes 
the diverging terms and thus prevents non-Markovianity being washed out in the higher 
temperature regimes. This is a feature specifically due to the spatial nature of the 
double well qubit, making this model well suited to realistic temperatures. 

In summary, the double well probe qubit model provides an ideal system for 
observing non-Markovianity in an atomic system, and for exploiting the Markovian 
to non-Mar kovian crossover to probe the BEC environment. In addition to involving 
systems that are straightforwardly combined experimentally, it is an example of a system 
showing a Markovian to non-Markovian crossover in accessible parameter ranges. We 
have shown here that in addition to the above advantages, this system is also robust to 
temperature effects, with the measured quantity maintaining its size for temperatures 
up to and, in some cases, beyond those necessary for experimental realization of these 
systems. 
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